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Abstract 

We calculate the leading quark mass corrections of order m q log(m q ), 

3 II 

m q and m q to the vector meson decay constants within Heavy Vector 
Meson Chiral Perturbation Theory. We discuss the issue of electromagnetic 
gauge invariance and the heavy mass expansion. Reasonably good fits to 
the observed decay constants are obtained. 



1 Introduction 



In this letter our aim is to determine within the Heavy Meson Effective Theory 
(HMET) the decay constants for the vector nonet. HMET was introduced |T|] as 
the non-relativistic limit of an interacting theory between vector mesons (heavy 
mesons) and a pseudoscalar meson background (light mesons). The theory is 
formulated in terms of operators involving the hadronic fields. The main reason 
to introduce such a formalism is to recover a well defined power counting in small 
masses and momenta. This is similar to the Heavy Baryon Chiral Perturbation 
Theory @. 

We list here the extra terms needed up to order p 3 for the vector decay con- 
stants. All the other relevant terms were already classified in |3[] and their coef- 
ficients estimated there. In principle these coefficients or Low-Energy-Constants 
(LEC's) should be determined directly from QCD. We simply fit their values to 
experiment, using Zweig's rule to limit the number of relevant constants, just 
as was done for the vector meson masses in ||. Here we fit 5 new data with 3 
new parameters. The poor experimental knowledge in the vector meson sector 
together with the rapid increase in the number of LECs when we introduce higher 
orders in the effective lagrangian is the main lack of this method. Some other 
recent papers using HMET for Vector mesons are ||, [5| . 

The vector decay constants are experimentally determined through the branch- 
ing ratios of p°, u, <fi — > e + e~ , fi + fi~ , via an electromagnetic current in the matrix 
element, or the branching ratios of t~ — > v T p~ , u T K*~ , via the vector part of the 
weak current. 

Corrections to the decay constants appear at order m q \og{m q ). We calculate 
here up to order m 3 ^ 2 . Two-loop contributions start appearing at order m 2 . 

We first discuss our notation and a few definitions. In Sect. |3| we discuss how 
gauge invariance can be used to connect terms at different orders in the heavy 
mass expansion. Using this we then proceed in Section |3] to list all the terms in 
the Lagrangian that are needed. Then we give our main results, the vector decay 
constants up to order p 3 in the HMET expansion. We then present the numerical 
results and compare with the experimental values. 

In the appendices we describe the slight extension needed for the weak currents 
and we quote approximate expressions for the vector isospin states. 

2 Definitions and Notation 

We define here our notation and the basis of chiral transformations. For a general 
introduction to Chiral Perturbation Theory see ||. 

Under SU(3)l X SU(3)r global chiral rotations the goldstone fields can be 
collected in an unitary matrix field, U (0) transforming as: 

U(4>) - g R U(cf>)gl g L xg R £ SU(3) L x SU(3) R (1) 
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For the chiral coset space SU(3)l x SU(3)r/ SU(3)y our choice of coordinates 
allows us to write: 



U (4>) = exp i 



« 2 (0), 



V2 



Jll _L JL 



7T^ 
_J[!L + JL 

K° 



K+ 
K° 



(2) 



where F 



92.4 MeV and u{<p) transforms as 



(3) 



with g) the so-called compensator field, an element of the conserved sub- 
group SU(3)v 

In the nonet case the annihilation modes for the effective vector meson fields are 
collected in a 3 x 3 matrix given by 



v/2 



K 



*o 



K 



*o 



(4) 



while its hermitian conjugate, W 7 ^, parametrizes the creation modes. In what 
follows we are only concerned with the transverse components of the fields, 
i.e. v ■ W = v ■ W> — 0, v is the chosen reference velocity for the heavy Vector 
meson. The 'longitudinal' component of can easily be expressed in terms of 
the transverse components. This is due to the fact that a massive vector meson 
has 4 components, but only 3 degrees of freedom. This should be understood in 
the remainder. 

Under chiral symmetry the effective vector fields transform as 



W„ 



h(<P,g)W,h\<P,g) Wl -> h(<P,g)Wlh\<P,g). 



(5) 



As has been mentioned already, our purpose is to compute the vector decay 
constants in the effective theory. They are defined through the matrix elements: 



(0 | qa^j I W) = F w e^ 



(6) 



for a vector meson state normalized to 1 and with momentum my/v^. mw is the 
mass of the relevant meson and e M its polarization vector. 

In order to introduce photons (we extend the formalism to weak currents in 
the App. [A]) we use the external field formalism ||, coupling the pseudoscalar 



fields to external hermitian matrix fields, a M , t> M defined as: 



eQA? + ... 
eQAf + ... 



(7) 
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where Q is the diagonal quark charge matrix, Q = (l/3)diag(2, —1, —1). 

This inclusion of a M , fields promotes the global chiral symmetry to a local 
one, allowing thus to define a covariant derivative and a connection: 

D ll U = d li U-irJJ + iUl lt , r M = i(« t [9 M -ir> + «[9 M -iZ>t). (8) 

Instead of using the r M and 1^ fields, we will use the combination 

Q± = e (u ] Qu ± uQu ] ) , (9) 

where Q± transforms as Q± — > h((f), g)Q±h((f), gy under chiral transformations. 

Notice that the insertion of external field (photon or a weak current) does not 
modify the power counting, i.e 0(e) = 1. 



3 Gauge and reparametrization invariance 

In this section we discuss briefly the constraints from reparametrization invari- 
ance and the additional constraints from gauge invariance on the HMET la- 
grangian. We will discuss it simply in terms of a single neutral vector meson and 
the photon. It is convenient to start from the relation between the relativistic 
and the effective fields [|], see also the discussion of 0: 

yix _ 1 ( e ~irn v vx W n + Jmvv-xyyiA + yylf Qq\ 

\f2m v \ ' v J \\ 

where the subindex v denotes the velocity of the heavy meson particle referred 
to an inertial observer, W v has momenta small compared to myv and satisfies 
v ■ W v = 0. The longitudinal component is suppressed by 1/my, see Sect. 4 
in ||. Instead, choosing a second reference frame related with the first by a 
Lorentz transformation we should get the same description of the physics. This 
fact relates the expression of the vector field in both frames: 

v -, w = v + q, W v v = e - imvq - x (W^, + w u q ■ W^j . (11) 

q is infinitesimal and satisfies v ■ q — since v 2 = w 2 = 1. This determines 
some 1/my coefficients in the chiral expansion, that can not be modified by 
non-perturbative corrections]/]. 

For the photon we need also to consider gauge invariance. The lagrangian has 
to be invariant under 

A^x) -> A^x) + d^e(x). (12) 

where e(x) is the gauge field. This transformation is in a particular frame, for 
fixed v. Relevant momenta are around ±myi> and 0, we therefore perform a 
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Fourier decomposition into low and high momenta in all the fields entering in the 
gauge transformation: 



A^x) = e- imvV - x A^x) + e imvV - x Al(x) + A^x) 
e( x ) = e- imvv - x e(x) + e imvv - x e\x)y + ~e(x) (13) 

where A M (x), A^x), e(x) and e(x) are all low momentum fields. As was men- 
tioned in || a similar decomposition allows to take into account properly the low 
momentum component for the vector meson field. 

Since the (low momentum) electromagnetic U(l) is a subgroup of (the low 
momentum) SU(3)y, the effect of e(x) and A^ will be ignored in what follows, it is 
treated by the covariant derivatives defined above. Contrary the high momentum 
of the electromagnetic field, U(l), is not included in that group. In order to 
obtain a gauge invariant lagrangian under this high momentum subgroup the 
electromagnetic field should transform as 

A^(x) -> A^(x) - im v v^e(x) + <9 M e(x). (14) 



The transformation in Eq. ([Tj) is the equivalent of Eq. (|il|) for the gauge 
invariance. As Eq. (11) it determines higher order coefficients in the 1/my 
expansion. 

For instance, if we take the following toy lagrangian: 

£ = Wl(A» + a 1 d»(v-A)+a 2 (v-d)A>^)+h.c. 

= W\v v (V > + a x d"A v + a 2 d v AA + h.c. (15) 

a high momentum gauge invariance transformation Eq. ( |T4"D implies: 

i i 

ai = , a 2 = • (16) 

my my 

those coefficients cannot by modified by non-perturbative corrections. Repa- 
rametrization invariance then requires in addition the (separately gauge invariant) 
term 

Xrd v W^ (» - & i A v \ + h.c. (17) 
m v v ' 

with the coefficient fixed. Notice that this is precisely the combination of terms 
that the relativistic term d v V" (d u A" — d"A u ) produces, taking into account v ■ 
W v = 0. So terms that are not gauge invariant at first sight can be made gauge 
invariant by adding terms of higher order in the heavy mass expansion. We will 
work in a fixed gauge, v ■ A = to avoid this complication. 
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4 The effective lagrangian 



To construct the relevant terms in the lagrangian we use both 

v-W v = andv-i = 0, (18) 

corresponding to the temporal gauge. Together with Eq. ([5]) which involves 
vertices with photons and pseudoscalars, we also use the following building blocks: 

X± = u*x u ^ =t u X^ u i u \i — wdpUu', (19) 

where x contains the quark mass matrix, x — 2I?odiag(m n , m^, m s ). 

We now construct the most general structure involving pseudoscalar, vector 
meson and photon fields which should be invariant under Lorentz transforma- 
tions, chiral transformations charge conjugation, parity and time reversal. 

Making use of the lowest order equation of motion, v ■ DW^ = 0, which 
eliminates terms removable by field redefinitions, and all other constraints, we 
have the following non anomalous lagrangian to leading order in 1/N C , with N c 
the number of colors and (B) = tr(£>): 

A = X 1 (WlQ + )A (i + A 2 <Wj{x +> + h.c. (20) 

And to next order in 1/N C : 

A = Aa(Wjg+>(x+>> + A 4 (Wt)( x+ Q + )> + h.c. (21) 

Notice that in the three flavour case, terms involving only a trace over Q± van- 
ishes, so they never appear. 

Terms at next-to-leading order in N c are Zweig rule suppressed and we will treat 
the terms in Eq. (|2lD as 0(p A ). 

At lowest order the odd intrinsic parity sector of the lagrangian is given by: 

£ 3 = w{{Wl W u }u a )v^ al3 + zA 5 ((Wl{Q+, u a })A v v^ + h.c.) . (22) 



In Eq. (p(l|), Eq. ( [21] ) and Eq. (|22|), all the coupling constants are real 
numbers^. 

In fact, reparametrization invariance requires the presence of higher order 
terms proportional to g and A5. They only contribute to the vector decay con- 
stants at order p A . The connection between a relativistic formulation and the 
present one can be done as in || |j but the external fields need to be split up as 
was done for the photon field in Sect. §. 



1 We need to use C, P and T separately here to prove this. The HMET is not a relativistic 
field theory so the CPT theorem is not valid. C connects and WT. It is only using both 
the requirement of a hermitian lagrangian and T, which also connects with W^, that we 
can conclude that the Ai are real. 
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1 ' V 

=>=0 > V 0r > < > (g) 



Figure 1: The three effective diagrams contributing to the width. A double line 
is a Vector Meson, a dashed line a pseudoscalar meson and the circled cross a 
vertex from £1,2,3- 

5 Calculation of the Vector Decay Constants 

The Vector Meson leptonic widths are given by 

SyraL F w 



T (W - / + r) = Hr^-f «, + 2mf)Jm w - 4m? (23) 



u m w 

for a |9j<7j| hadronic system which decays via W — > 7 — > / + / _ with a eF^A^W^ 
effective coupling. Otherwise if a weak decay current is involved (i.e. a r lepton 
decay), the width is determined by 

T( r _> Wu) = \ Vckm \ 2 ° 2 f F w K - "4) 2 K + 

87r m w m\ 

where Gf = V2e 2 / (8 sin 2 QwMyy) is the Fermi constant and an effective coupling 
F w [e / (2^/2 sin 9 w)Vc km] has been used. Vckm is equal to Kd, V us for 

r + decay to p + v T , K* + v T respectively. rri\y is the mass of the vector meson, 
Mw = 80.3 GeV is the mass of the weak vector boson. The only difference with 
the definition is that we use the electromagnetic current for the neutral vector 
bosons. 

To determine inside the effective theory the Fw constants of Eq. (p3|) and ([24]) 
one has to compute the diagrams depicted in Fig. [l] and in addition the contribu- 
tion coming from the Vector Meson wave function renormalization (w.f.r.). For 
that purpose we first define the physical vector meson basis, to be used in what 
follows (T ext ), fully diagonalizing the two point Green function (G 2 (p 2 , m 2 hys )) at 
the 1-loop level, this contribution is found using the results of Ref. || for G 2 up 
to 0(p 4 ). This is sufficient to compute the w.f.r. as defined by 

Z v = — G 2 (k w , 

iriphys) \on— shelh 

(25) 

to order p 3 , they include the electromagnetic corrections of [10|. Here kw^ = 



Pwv ~ m v v iJ. an d Pw = m wi kw is the momentum in the HMET. 

The other contributions are found by direct calculation using Eq. (p0|) and 
Eq. (|H). We define 

X°=X+L=o and Q° + = Q + \ U=0 (26) 
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Figure 2: One of the relativistic diagram contributing to the decay constant. 

to be the quantities Q + and x with u = 0. The extension of Q + needed for the 
charged case is defined in (|A.4j) . The contribution of the tree diagram, the first 
diagram of Fig. |], to the decay constant corresponding to T£ xt is 

MT e a xt Q° + ) + A 2 <T^{ X ° , Q° + }> . (27) 
The tadpole type (second) diagram is given by: 

E ^(TUT M ,[T M \Q + ]})m, (28) 

M=l,8 T 

and finally the sunrise type diagram (last diagram): 

E E ^({^, ^ t }r M > (^ t {^, r Mt M(Pe, tl Am C) m M ) , (29) 

M=l,8c=l,9 t 

where the T int basis was defined in M and the fiM , K(p ext , Am c , tum) functions 



are defined in Eq. ( |B.5|) of App. (0). The full result is then given by the 



sum of equations (pTj ), (|28|) and (P9|) divided by the square root of the relevant 
Zy. We have kept also some higher order terms required by reparametrization 
invariance. This corresponds to using instead of the function K(p ext , Am c , tum) 
the full combination [G, fi] defined in ( |B.6| ). 

We have checked that the non-analytical pieces of the relativistic diagram in 
Fig. |^ are fully recovered by the second diagram in Fig. |T] as was explicitly shown 
in 01 for the scalar form-factor case. 



6 Numerical results and conclusions 



From the decay widths of Eqs. (|23|j2^ ) we obtain the experimental results of 



the second column of Table |l[ The Cabibbo-Kobayashi-Maskawa mixing angles, 



lifetimes, masses and branching ratios are taken from the particle data book [ |TT 
The naive prediction, using all the pure isospin states, is 

\/2F p o = 3V2F W = -3F = F p+ = F K *+ . (30) 



2 We need to extract a factor of e, respectively ej '(2y/2 sin 6 W)^C KM, compared to the defi- 
nition of Q° + in (p6l) to obtain the decay constant Fw 
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Scenario 


exp 


LO 


III 


IV 


VI 


lOOAi (GeV 3 / 2 ) 




6.42 


6.97(36) 


7.55(39) 


7.74(38) 


100A 2 (GeV- 1 / 2 ) 




0.00 


-0.26(50) 


0.40(49) 


1.21(69) 


100A 5 (GeV 1 / 2 ) 




0.00 


3.56(75) 


5.12(71) 


7.24(1.43) 


F p0 (GeV 3 / 2 ) 


0.0946(23) 


0.0908 


0.0951 


0.0955 


0.0937 


F u (GeV 3 / 2 ) 


0.0287(5) 


0.0303 


0.0286 


0.0278 


0.0282 


-F+ (GeV 3 / 2 ) 


0.0565(10) 


0.0428 


0.0563 


0.0544 


0.0534 


F p+ (GeV 3 / 2 ) 


0.1284(10) 


input 


0.1275 


0.1234 


0.1280 


F K *+ (GeV 3 / 2 ) 


0.1447(45) 


0.1284 


0.1460 


0.1571 


0.1554 



Table 1: The experimental Decay Constants and various good fits. The column 
labeled LO (lowest order) is Eq. ([50) with F p + as input. For an explanation of 
the scenarios and the values of the other constants see Table 1 in J|. 

This is satisfied to about 5% for p + ,p°,uj, to about 13% for p + ,K* + and about 
32% for p + , (p. The signs we have fixed to agree with (p0|). 

The input parameters used are scenarios III, IV and VI from Ref. 0. Ill 
and VI were fits to the masses only for two different values of g, one high and 
one low, fit IV also included the p — uo mixing in the input but otherwise as fit 
III. We find a reasonable fit to all the decay constants for reasonable values of 
Ai 2,5. The higher order corrections are also reasonable, below 40%. The fit for 
scenario IV is worse for the following reason: Using the mixings from || for the 
u, including p 4 effects, the contributions at tree level from the Ai term essentially 
cancel, leaving the loop diagrams of Fig. [I] and A 2 as the main contributions. 
This makes the predictions for the us somewhat unstable. We also have a rather 
large w.f.r. factor for the <fi. The total size of the higher order corrections can 
be judged by comparing the results from Eq. (|3~0"D, column LO, with those of the 
three scenarios. 

What we have minimized in order to get the values of Ai^.s in Table [l] is 

Ffit = Ew=i,5 (\Fw eo /Fw P \ ~ l) • The error on the Aj corresponds to changes in 
Ff it by about 0.01, i.e. at most 10% for an individual F w , minimizing the other 
two Aj at the same time. 

In conclusion we have calculated the corrections to the vector decay constants 
in heavy vector meson chiral perturbation theory and found acceptable fits to all 
the measured ones. We have determined 5 observables in terms of 3 parameters. 

A Weak currents 

In this appendix we give the main features to incorporate charged weak current 
effects to our formalism. The neutral weak current is not phenomenologically 
relevant at present. To do so one needs to extend the left current defined in Eq. 
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©to: 



Ifi — cQA^ + 



(a.i; 



a/2 sin 6*iy 

where sm(9w) is Weinberg's angle, W M parametrizes the spin-1 gauge boson 
fields - it creates an W + gauge boson field and destroys an W~ one - and we 
have introduced the T matrices in terms of the relevant Cabibbo- Kobayashi- 
Maskawa factors 



( o v ud v us 



V o o o 



7i 



4 





v ud 

V n , 



(A.2) 



As in the QED sector, the requirement of not breaking chiral symmetry force us 
to split the field in different components in the momenta space according to 



VV^ + e^^VVL. (A.3) 



The inclusion of the charged currents is now achieved by replacing 



eQ+A, 



by 



eQ+A^ + 



a/2 sin 9 W 



u{T+W. 



(A.4) 



in lagrangians Eq. (|20|) , Eq. (]2l|) and Eq. (0). Where now P violation is 
allowed. 



B Vector Decay Constant Contributions 

In this appendix we show the formulae for the contribution coming from the 
effective diagrams of Fig. p]. Where we take the approximation of non-diagonal 
fields inside T ext and T int , i.e. we use the pure isospin 1 state for p°, pure isospin 
for uj and and the as the pure strange vector state, i.e. we neglect here 
p° — uj — mixing. 

For the tree level contribution we find 

8 v^2 

F p o = v^Ai + —^-B \ 2 (m d + 2m u ) 

Fu, = - ^-^B X 2 (m d - 2m u ) 

F<j> = — -Ai — —B m s X 2 

F p + = 2Ai + 8B (m u + m d )X 2 
F K *+ = 2\ 1 + 8B (m u + m s )\ 2 . (B.l) 
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For the tadpole type diagrams, we find 
Ai 



F p0 



ll K + 

PK+ 



F + = — 
p F 2 \ 2 



Ai 
V2F* 

Fl 



H h ^71 + + C /U^O + S fijj 



Ai 



+ + ^ + (v 3c + s) -i + (c - V3s) — . 



2 fin 



(B.2) 



And finally the contribution coming from the sunrise diagram 



F p0 



A 



( cs 



V3V6 y/2 



( C 



cs 



3V2 3VQ 



hP } 



1 r „j V^r^O r^*0l 



+ 



/ C 1 



CS 



\V2 3VQ 



3V2 3VQ 



( / c cs \ / cs s \ 



1 V^r^Q r^*0l 



c,s 



V9V2 



cs 



V3\/2 >/6 



3 I 3 



a{|W] + \iK + ,K*°] + \[K\K* + } + [tt+ u;] + j[^°,P + ]} 



with 



A = ^^ 5 , 0086* = c and sin# = s. 



(B-3) 



(B.4) 



In addition to Eq. ( |B.1| ), Eq. ( |B.2| ) and Eq. ( |B.3| ) one has the w.f.r. terms. 
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We have defined the following integrals 

d d q q^q u 1 



Iflr, 



(2ir) d v ■ q — 00 + irj q 2 — m 2 + it] 
d d q 1 im 2 



(27r) d q 2 -m 2 16tt 2 
with A = 1/e — 7 + log(47r) + 1 and d = 4 - 2e, and 



(B.5) 



[G, O] = ( 1 + ^ - ) K(m G , mn - mv -p. v ) (B.6) 



Here 77 and ttq stand for the physical pseudoscalar fields. 
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